Multiplexity-facilitated cascades in networks 
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Elements of networks interact in many ways, so modeling them with graphs requires multiple 
types of edges (or network layers) . Here we show that such multiplex networks are generically more 
vulnerable to global cascades than simplex networks. We generalize the threshold cascade model 
[D. J. Watts, Proc. Natl. Acad. Sci. U.S.A. 99, 5766 (2002)] to multiplex networks, in which a node 
activates if a sufficiently large fraction of neighbors in any layer are active. We show that both 
combining layers (i.e., realizing other interactions play a role) and splitting a network into layers 
(i.e., recognizing distinct kinds of interactions) facilitate cascades. Notably, layers unsusceptible to 
global cascades can cooperatively achieve them if coupled. On one hand, this suggests fundamental 
limitations on predicting cascades without full knowledge of a system's multiplexity; on the other 
hand, it offers feasible means to control cascades by introducing or removing sparse layers in an 
existing network. 

PACS numbers: 89.75.Hc, 87.23.Ge 



When choosing which products to buy, ideas to adopt 
and movements to join, people are influenced by friends, 
colleagues, family, and other types of contacts. Such 
influence along multiple channels is frequently non- 
additive: Just one type of relationship often suffices to 
convince someone to change behavior [1 . Banks also 
interact in many ways — through balance sheet claims, 
derivatives contracts and reliance on credit lines — which 
collectively and nonlinear ly cause cascades [2 . At 
a broader scale, countries interact not only through 
trade [3] but also through investment and lending [4]. 
Models with just one type of edge [5 cannot capture 
this non-additive influence along multiple channels. In- 
stead, one needs graphs that explicitly contain multiple 
types of edges (or network layers), called multiplex net- 
works [11 [6l [7] . Such multiplex networks provide a com- 
plementary framework to the growing body of works on 
interacting and interdependent coupled network systems 

Here we study the impact of such network multiplexity 
on cascade dynamics in the threshold model introduced 
by Watts [13 . In this stylized model of, for example, con- 
tagious behavioral adoption in a social network [Ml [15] , 
people join the growing movement if a sufficiently large 
fraction of their friends have. Similarly, banks default if 
sufficiently many debtor banks default [16]. Specifically, 
nodes exist in one of two states, active and inactive. Each 
node independently draws a (frozen) threshold r G [0, 1] 
from a probability distribution Q{r). A node of degree 
k activates if its fraction m/k of active neighbors ex- 
ceeds its threshold r. Of particular interest are so-called 
global cascades, in which a finite fraction of the infinite 
network becomes activated from a vanishingly small frac- 
tion of initially active seeds. A key lesson from previous 
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studies is that, for a given distribution of thresholds, net- 
work connectivity constrains global cascades [13, 17 . If 
it is too sparse, a network lacks a giant component and 
the connectivity needed for a global cascade; if it is too 
dense, a network likely cannot surround nodes with suf- 
ficiently many active neighbors. Various generalizations 
have since been introduced and studied p!6ti2Q] . 

In this Rapid Communication we generalize Watts' 
threshold model [13] to multiplex networks, in which 
nodes activate if a sufficiently large fraction of neighbors 
in any layer are active. To motivate this formulation, 
note that in many situations what matters is the influ- 
ence from one layer alone. For example, a large fraction 
of colleagues recommending a certain smartphone appli- 
cation may convince someone to use it. Similarly, the 
default of sufficiently many loans may suffice to depress 
a country's trade, and vice versa. To be specific, in a 
network with two layers (a duplex network), a node with 
ki and k2 many neighbors in layers 1 and 2, respectively, 
with mi and m2 of those neighbors active, itself activates 
if mi/ki or m2//^2 exceeds the node's threshold r. We 
denoted this multiplex model the 102 model. For com- 
parison, we also consider the simplex network that has 
the same topology but that ignores multiplexity, denoted 
the 1©2 model (Fig.[l]). 

The central result of this Rapid Communication is 
the greater ease of cascades in multiplex networks. We 
demonstrate the effect of multiplexity on cascades in two 
scenarios. First, given a singe-layer network [denoted 
Layer 1 in Fig.flla)], one might realize that another kind 
of interaction (Layer 2) plays a role. In this case, we 
combine a second layer (Layer 2) with the existing one 
to form the duplex network l(g)2 [Fig. [ijb)]. Alterna- 
tively, one might realize that a given network in fact con- 
sists of multiple channels of non-additive interactions, so 
the simplex network [Fig. [TJc)] is split into two layers 
to form the l(g)2 network [Fig.jljb)]. We find that both 
combining layers and splitting into layers facilitate global 
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FIG. 1: (Color online) Construction of the multiplex 102 net- 
work (b) by combining two layers (a), the solid green (top) 
and dashed red (bottom) interactions. The simplex 102 net- 
work (c) ignores the types of interaction (solid gray edges). 



cascades. Layers that in isolation have too much or too 
little connectivity to achieve global cascades can cooper- 
atively achieve them if they are multiplex-coupled. We 
demonstrate this analytically and using simulations, and 
we conclude with generalizations to networks with three 
or more layers. 

We begin by extending the theory of Ref. [17 to a 
duplex network with layers of locally tree-like random 
graphs on the same set of N nodes. Every node has two 
independent degrees, ki and /c2, equal to its numbers 
of neighbors in layers 1 and 2, respectively. The mean 
fraction of active nodes in the stationary state, called 
the mean cascade size p, equals the probability that a 
randomly chosen node is active. This probability can be 
obtained by approximating the network as a tree with the 
chosen node as its root and by considering the cascade 
of activations toward the root [17 . Given the initial seed 
fraction po^ P foi" a duplex 10)2 network is given by 

p = po + (l-po) Pk^Pk2^ 

/C1+/C2>1 

X iz E <\('/^^)st('/S^)^^v'4, (1) 

mi=0 7712=0 

where {q^ , ) is the fixed point of the recursion 
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begun from q^^ = po^ where, for distinct z, j G {1,2}, 

9'H<i'r^\e') = + (1 - po) E E 

X E E B^r,l-\q'^^)B%{qliy)F^^^X. (3) 
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Here p)^^ is the degree distribution of layer i; B^{q) = 
(^)q'^{l — q)^~^ is shorthand for the binomial distribu- 
tion; qn^ is the probability that a node n steps above the 




FIG. 2: (Color online) Theoretical cascade boundary given 
by the first-order cascade condition ^ (dashed green line) 
and numerically simulated mean cascade size p for po = 10~^ 
(color coded) on a duplex network oi N — 10^ with ER layers 
of equal mean degrees z. Inset: p \s z for threshold = 0.18 
and different po = 5 x 10"^ (red A), 10"^ (blue O), 5 x 10"^ 
(□), obtained from simulations (symbols) and from Eq. ([T]) 
(lines). 



leaves of the tree is activated by its children in the tree, 
conditioned on its parent in layer i being inactive; F^^^^^ 
is the response function, the chance that a node with ki 
neighbors in layer i (m^ of which are active) becomes 

active; and the factor kip)^^, / zi in ([3| is the probability 
that a degree- /c^ node lies at the end of a randomly cho- 
sen edge in layer z, where Zi is the mean degree in layer 
i. The response function for the multiplex 102 model is 



mi,m2 



if max(mi//ci, m2//c2) ^ ^, 
if max(mi//ci, m2//c2) > ^- 



(4) 







In this work all nodes have the same threshold R [i.e., 
Q{r) = S{r - R)]. 

To test the validity of the theory, we calculated p from 
jj) for duplex Erdos-Renyi (ER) networks [21] of layers 
with the same mean degree z. The calculated p as a func- 
tion of z with different seed sizes po are found to agree 
well with numerical simulations (Fig. [2j inset). As in 
the single-layer case [131 llZ] ? foi" given i?, global cascades 
occur for an interval of z between two transitions: a con- 
tinuous transition at small z for the emergence of global 
cascades (following the emergence of the giant connected 
component), and a discontinuous transition for the dis- 
appearance of global cascades, the location of which in- 
creases with pq. The plot of p in (i?, z)-parameter space 
(Fig. [2j main plot) shows the cascade region, the param- 
eters for which global cascades occur. 

The linear stability of the fixed point (<7^^\<7^^^) = 
(0, 0) of the recursion ([2| as po ^ gives a sufficient con- 
dition for global cascades, leading to the so-called first- 
order cascade condition [13l [ITl |19l [20] . In the multiplex 
102 case, this condition is that the maximum eigenvalue 
of the Jacobian matrix J of the recursion ([2| at the origin 
as Po ^ exceeds 1: 



Amax(J) > 1- 



(5) 
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The 2x2 matrix J is given by, from (3) and (4), 



L1/-RJ 

/Ci=l 



ki{ki 



Dpi' 



[l/Rj 



Jl2 = ^^Pk: 



(2) 



(6) 



and similarly for J21 and J22, where [-J denotes the floor 
function. As shown in Fig. |2] (main plot), this cascade 
condition ([5| closely approximates the boundary of the 
cascade region from simulations, providing a useful ap- 
proximation to the actual cascade region. Deviations 
from simulations occur because ([5| ignores activations 
by multiple active neighbors. (Including second-order ac- 
tivations improves the agreement [17]). For the exact 
boundary, one must use Eqs. ([T])-(|3|, as in Fig. |2] (inset). 

The following interpretation of the cascade condition 
elucidates how multiplexity facilitates cascades. The ma- 
trix J can be identified with the mean reproduction ma- 
trix of a two- type branching process [22] , with Jij repre- 
senting the mean number of offspring of type-j branch- 
ing from a node born through a type-z branching. Here 
a type-z branching corresponds to activation along the 
layer z by a single active node. The cascade condition 
is the supercriticality condition for this two-type branch- 
ing process. This branching process approximates the 
progression of actual cascades by ignoring activations by 
multiple active neighbors. Cascades in a single-layer net- 
work are approximated by a single-type branching pro- 
cess, so the cascade condition becomes Jii{R, z) > 1 [T3] . 
For the l(g)2 network, the cascade condition reads Amax = 



> 1. Since 



( Jll + J22) + V ( Jll - J22)^ + 4Ji2 J21 

this Amax > niax(Jii, J22), the l(g)2 cascade region con- 
tains the cascade regions of either of its layers in isola- 
tion, for any degree distributions. Thus what enlarges 
the cascade regions is the presence of an additional acti- 
vation channel, represented by the off-diagonals Jij. The 
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FIG. 4: (Color online) Cascade regions with = 0.18 for (a) 
the union of two independent ER layers with mean degrees zi 
and Z2 (denoted 1U2), (b) the simplex (102) ER network with 
mean degree zi + ^2, and (c) the duplex (102) network with 
ER layers of mean degrees zi and Z2. The dotted green lines 
are the cascade boundary obtained from Eq. (jsj); numerically- 
simulated p (po — 10~^) is color-coded. Also shown are addi- 
tional cascade regions in the multiplex cases for (d) combining 
and (e) splitting. 
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FIG. 3: (Color online) Both combining layers (a) and split- 
ting into layers (b) facilitate cascades, (a) Cascade regions 
from cascade condition ([5]) for single-layer (yellow, brighter) 
and duplex (red, darker) network with ER layers of equal 
mean degree z. (b) Same plots for a simplex ER network of 
mean degree 2z (yellow, brighter) and corresponding duplex 
network (red, darker). Overlapping regions appear orange 
(medium brightness). 



two types of activation channels promote each other's 
activations in a cooperative, positive- feedback manner, 
collectively facilitating cascades in multiplex networks. 

As an explicit example, we consider duplex networks 
with ER layers of equal mean degrees under the two sce- 
narios of multiplexity discussed above. First, a second 
ER layer with the same mean degree z is combined with 
an existing ER layer with mean degree z to form the 102 
network (of total mean degree 2z), akin to considering 
lending as well as trading relationships among countries 
in the global economic system. Second, an ER graph 
with mean degree 2z is randomly split into two ER lay- 
ers with equal mean degrees z to form the 102 network, 
akin to distinguishing social influence among colleagues 
and among friends in a social network. Explicit evalua- 
tion of the cascade condition (5) yields, in both cases, an 
enlarged cascade region for the 102 network (Fig. |3|. 

A multiplex network with statistically distinct lay- 
ers supports even more nontrivial cascades. Both com- 
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FIG. 5: Cascade condition of ^-plex networks of ER layers 
each with mean degree z, for i = 1, 2, 3, 4 and 10. Note that 
the cascade region extends to > 1/2 for > 4. We show 
white, dashed hnes on boundaries of regions for visual clarity. 

bining two ER layers with mean degrees zi and Z2 
[Figs. [4]^ a), |4]^c), and |4jd)] and splitting an ER net- 
work into two layers with mean degrees Zi and Z2 
[Figs. [4]^b), [4]^c), and [4]^e)] enlarges the cascade region. 
The additional cascade regions in the multiplex case 
[Figs, ^d) and |4je)] highlight the cooperative effect of 
multiplexity. For these parameters, each layer is too 
sparse or too dense to achieve global cascades, but they 
cooperatively achieve them when multiplex-coupled. Of 
particular interest is when one of 2^1, ^2 is too small (< 1), 
the other too large 6) to support global cascades in 
isolation. This presents a way to control cascades in a 
system by introducing or removing sparse layers below 
the percolation threshold, which may be more feasible to 
implement than perturbing the existing, dense network. 

We verified a similar enlargement of cascade regions 
for networks with layers of broad degree distributions by 
using the static model of scale-free graphs, which allows 
a variable mean degree [23]. Furthermore, we checked 
that short loops introduced by multiplexity can be ne- 
glected [19 for the large, sparse networks considered 
here. 

Finally we generalize to multiplex networks with i > 2 
layers. Extending Eqs. ([l])-(|6| to i layers is straight- 
forward. Combining more layers further facilitates cas- 



cades for larger z and R (Fig.[5|. Notably, introducing a 
fourth layer permits global cascades even for thresholds 
R > 1/2, which Morris [2j proved cannot occur in sim- 
plex networks. This can be understood from the cascade 
condition as follows. For R > 1/2, only activations of 
degree- 1 nodes contribute to the cascade condition (5) 
and (6). For multiplex networks with i ER layers with 
equal mean degree the first-order cascade condition 
for R > 1/2 thus becomes {£ — l)ze~^ > 1, which can 
be satisfied by a non-empty interval of z for ^ > 4. This 
suggests that even people difficult to persuade {R> 1/2) 
to buy a new device, for example, may all buy one if they 
participate a little ^ 1) in many social spheres > 4). 

To conclude, the interplay among multiple kinds of 
interactions — the multiplexity — can generically increase 
a network's vulnerability to global cascades in a thresh- 
old model. Interestingly, layers unsusceptible to cascades 
can cooperatively become susceptible when coupled. The 
impact of multiplexity on network dynamics is expected 
to be widespread [SHIS [25] . other binary-state, mono- 
tonic, threshold cascade models such as bootstrap perco- 
lation [26], a similar analysis can be readily applied, while 
its impact on more complicated dynamics |27| remains 
largely unexplored. Our results suggest a double-faceted 
picture for cascade prediction and control. On one hand, 
one faces fundamental limitations on predicting cascades 
without full knowledge of a system's multiplexity; on the 
other hand, multiplexity offers a feasible tactic to enable 
or hinder cascades by introducing or removing sparse lay- 
ers, respectively. For instance, advertising may become 
more effective with every new medium, while banks may 
grow more vulnerable with every new lending mechanism. 
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